Dark matter annihilation: the origin of cosmic gamma-ray background at 1—20 MeV 



Kyungjin Ahn and Eiichiro Komatsu 
Department of Astronomy, University of Texas at Austin 
1 University Station, C1400, Austin, TX 78712 



The origin of the cosmic 7-ray background at 1-20 MeV remains a mystery. We show that 7-ray 
emission accompanying annihilation of 20 MeV dark matter particles explains most of the observed 
signal. Our model satisfies all of the current observational constraints, and naturally provides the 
origin of "missing" 7-ray background at 1-20 MeV and 511 keV line emission from the Galactic 
center. We conclude that 7-ray observations support the existence of 20 MeV dark matter particles. 
Improved measurements of the 7-ray background in this energy band undoubtedly test our proposal. 
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PACS numbers: 95.35.+d, 95.85.Nv, 95.85.Pv 



What is the origin of the cosmic 7-ray background? 
It is usually understood that the cosmic 7-ray back- 
ground is a superposition of unresolved astronomical 7- 
ray sources distributed in the universe. Active Galac- 
tic Nuclei (AGNs) alone explain most of the background 
light in two energy regions: ordinary (but obscured by 
intervening hydrogen gas) AGNs account for the low- 
energy (< 0.5 MeV) spectrum [ESQ, whereas beamed 
AGNs (known as Blazars) account for the high-energy 
(> 20 MeV) spectrum 0,0, El- There is, however, a gap 
between these two regions. While historically supernovae 
have been a leading candidate for the background up to 

4 MeVUHUGH, recent studies[H|l2 show that the su " 
pernova contribution is an order of magnitude lower than 

observed. The spectrum at 4-20 MeV also remains unex- 
plained (for a review on this subject, see [1^1. It is not 
very easy to explain such high-energy background light 
by astronomical sources without AGNs or supernovae. 

So, what is the origin of the cosmic 7-ray background 
at 0.5-20 MeV? On energetics, a decay or annihilation of 
particles having mass in the range of 0.5 MeV < mx 
20 MeV would produce the background light in the de- 
sired energy band. Since both lower- and higher-energy 
spectra are already accounted for by AGNs almost en- 
tirely, too lighter or too heavier (e.g., neutralinos) par- 
ticles should be excluded. Is there any evidence or rea- 
son that such particles should exist? The most com- 
pelling evidence comes from 511 keV line emission from 
the central part of our Galaxy, which has been detected 
and mapped by the SPI spectrometer on the INTEr- 
national Gamma-Ray Astrophysics Laboratory (INTE- 
GRAL) satellite [3, This line should be produced 
by annihilation of electron-positron pairs, and one of the 
possible origins is the dark matter particles annihilat- 
ing into electron-positron pairs fl6j. This proposal ex- 
plains the measured injection rate of positrons as well 
as morphology of the signal extended over the bulge re- 
gion. Intriguingly, popular astronomical sources such as 
supernovae again seem to fail to satisfy the observational 
constraints [13I . Motivated by this idea, in the previous 
paper 18] we have calculated the 7-ray background of red- 
shifted 511 keV lines from extragalactic halos distributed 
over a large redshift range. We have shown that the an- 



nihilation signal makes a substantial contribution to the 
low-energy spectrum at < 0.511 MeV, which constrains 
mx to be heavier than 20 MeV in order for the sum of 
the AGN and annihilation contributions not to exceed 
the observed signal. 

In this paper, we extend our previous analysis to in- 
clude continuum emission accompanying annihilation. 
The emerging continuum spectrum should of course de- 
pend on the precise nature of dark matter particles, 
which is yet to be determined. Recently, an interest- 
ing proposal was made by |l9j |: radiative corrections 
to annihilation, XX — > e + e~, should lead to emission 
of 7-rays via the internal bremsstrahlung, the emission 
of extra final-state photons during a reaction, XX — > 
e + e~7. They have calculated the spectrum of the in- 
ternal bremsstrahlung expected for annihilation in the 
Galactic center, compared to the Galactic 7-ray data, 
and obtained a constraint on mass as mx ^ 20 MeV. A 
crucial assumption in their analysis is that the cross sec- 
tion of internal bremsstrahlung is linearly proportional 
to the annihilation cross section, and the constant of 
proportionality is independent of the nature of annihi- 
lation, as is found for related processes [2(1 0, |2^, • 
More specifically, they assumed that the cross section of 
XX — > e + e~7 would be calculated by that of e + e~ — > 
/j + /i~7 with the muon mass replaced by the electron 
mass. Although the equivalence between these two 
processes/cross-sections has not been demonstrated as 
yet, we adopt their procedure into our calculations. 

We calculate the background intensity, l v , as |24| 



l _ c_ l dzP u ([l + z]v,z) 
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H(z)(l + z) 4 



(1) 



where v is an observed frequency, H(z) is the expansion 
rate at redshift z, and P v {y,z) is the volume emissivity 
(in units of energy per unit time, unit frequency and unit 
proper volume): 



P v = \^v(av)n\ 



4a g(y) 



(2) 



where a ~ 1/137 is the fine structure constant, nx 
is the number density of dark matter particles, and 
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(av) is the thermally averaged annihilation cross sec- 
tion. To fully account for WMAP's determination of 
mass density of dark matter pfij. flxh 2 — 0.113, by 
cold relics from the early universe, one finds (av) = 
[3.9,2.7,3.2] 1(T 26 cm 3 s" 1 for m x = [1, 10, 100] MeV, 
respectively (e.g., see Eq. [1] in ,26]). We have assumed 
that (av) is velocity- independent (S-wave annihilation). 
One might add a velocity-dependent term (such as P- 
wave annihilation) to the cross-section; however, such 
terms add more degrees of freedom to the model, mak- 
ing the model less predictable. While Bcehm et al.flfij 
argue that the S-wave cross section overpredicts the 7- 
ray flux from the Galactic center, we have shown in the 
previous paper0] that it is still consistent with the data 
for mx ^ 20 MeV and the Galactic density profile of 
p oc r~ 0A or shallower. (We shall discuss an issue re- 
garding the density profile later.) Finally, a dimension- 
less spectral function, g(v), is defined by 
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where s' = Amx(mx — his). This function is approxi- 
mately constant for hv < mx , and then sharply cuts off 
at hv ~ mx- Thus, one may approximate it as 
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FIG. 1: Cosmic 7-ray background from dark matter annihi- 
lation. The dotted, short-dashed, and long-dashed lines show 
mx = 10, 20, and 50 MeV, respectively. The curves which 
sharply cut off at 511 keV represent background light from 
line emission 18], while the others which extend to higher en- 
ergy represent the internal bremsstrahlung. 



for the sake of an order-of-magnitude estimation. (Note 
that we have also assumed mx 3> m e .) 

Since the number density is usually unknown, we use 
the mass density, px = nx/mx, instead. After multi- 
plying by v, one obtains 
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where {p x ) is the average of p x over proper volume 
at z, and Cx(z) = (px) z / (Px) z ^ s the dark matter 
clumping factor. (We have used (px) z — 10.54 Qxh 2 (l + 
z) 3 keV cm -3 .) While equation (JjjJ is exact, one may ob- 
tain a better analytical insight of this equation by using 
the approximation to g(v) (Eq. [I]), 
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where y = hv(l + z)/mx- Here, we have also assumed 
that the integral is dominated by 1+z ^> (Cl^/Cl m ) = 2.3. 

We follow the method developed in our previous paper 
[l8j for calculating the clumping factor of dark matter, 
Cx(z). We have shown that C x (z) at z < 20 is approx- 
imately a power law, 



Cx(z) = C x (0)(l + z)- 



(7) 



and depends on adopted dark matter halo profiles. 
For example, a cuspy pro file such as the Navarro-Frenk- 
White (NFW) profile [13, px[r) oc r~\ gives C x (0) ~ 
10 5 and ~ 1.8, while a flat profile such as the Trun- 
cated Isothermal Sphere (TIS) [2^, Px(r) oc r°, gives 
C x (0) ~ 10 3 and ~ (see Figure 2 of Q). Using a 
power-law evolution of Cx(z), one obtains an approxi- 
mate shape of the spectrum as 



hv\n(2mx/m e ) 



( hv 
\m x 



/3-3/2' 



i'J r x ; ; — — — t— 1 — I ] d(mx — hv), 

(8) 

for m x > m e . If < 3/2 (e.g., TIS), vl v oc 

(hv)P~ 1 / 2 (himx)/m x +1 ^ 2 i9(mx — hv), whereas if > 
3/2 (e.g., NFW), vl v oc hv[ln(2m x /m e )]/m 2 x 'd(mx - 
hv). Note that the shape of the spectrum becomes in- 
sensitive to halo profiles for the latter case (while the 
amplitude still depends on profiles). 

Henceforth we shall adopt the NFW profile as the fidu- 
cial model, as it fits the mean central halo profiles in nu- 
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FIG. 2: The total cosmic 7-ray background produced by dark 
matter annihilation, AGNsQ, and Type la supernovae|l2jl. 
The dotted, short-dashed, and long-dashed lines show mx = 
10, 20, and 50 MeV, respectively. The supernova contribution 
depends on the observed supernova rate, and we consider the 
best-fit rate (upper panel) as well as the 2a upper limit (lower 
panel). The datapoints of HEAO-1 A4 MED^, SMM^, 
and COMPTEL^J experiments are also shown. 
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FIG. 3: The best-fit model of the cosmic 7-ray background. 
The model assumes (a) mx = 20 MeV, (b) the mean dark 
matter clumping factor is twice as large as predicted by the 
NFW profile (due to either a steeper profile or the presence 
of substructures), and (c) line emission is solely produced via 
positronium formation. The dashed lines show each contribu- 
tion separately. 



merical simulations well. Following the previous paper, 
we take into account a scatter in halo profiles by integrat- 
ing over a probability distribution of halo concentration; 
thus, our model effectively incorporates significantly less 
concentrated (such as our Galaxy) or more concentrated 
profiles than the average NFW. One might argue that 
our model based on the NFW profile is unable to explain 
7-ray emission from the Galactic center, which requires 
p oc r -0 ' 4 (or shallower). If desired, one might use this 
profile and recalculate the 7-ray background spectrum; 
however, we continue to use the NFW profile, assum- 
ing that our Galaxy is not a "typical" halo in the uni- 
verse. If there are so many more galaxies which obey 
the NFW profile, then the signal should be dominated 
by those typical halos. Of course, real universe does not 
have to be the same as numerical simulations, and one 
way to incorporate the uncertainty of halo profiles into 
our analysis would be to treat Cx(0) and (3 as free pa- 
rameters. We shall come back to this point at the end 
of this paper. Figure ^ shows the predicted cosmic 7- 
ray background from dark matter annihilation, including 
line|l8j| and continuum emission, for mx = 10, 20, and 
50 MeV. The shape of the internal bremsstrahlung is de- 
scribed well by the approximate formula (Eq. with 
/3 = 1.8. As expected, the continuum spectrum extends 
up to hv ~ mx, whereas line emission contributes only 
at < 0.511 MeV. 



Now let us add extra contributions from known astro- 
nomical sources and compare the total predicted spec- 
trum with the observational data. Figure compares 
the sum of dark matter annihilation, AGNs 3] and Type 
la su pern ovae[l2j with the data points of HEAO-1 [23, 
SMM|13, and COMPTELjp experiments. We find that 
mx ~ 20 MeV fits the low-energy spectrum^] and ex- 
plains about a half of the spectrum at 1-20 MeV. There- 
fore, the internal bremsstrahlung from dark matter an- 
nihilation is a very attractive source of the cosmic 7-ray 
background in this energy region. It is remarkable that 
such a simple model provides adequate explanations to 
two completely different problems: 511 keV line emis- 
sion from the Galactic center and missing 7-ray light 
at 1-20 MeV. (The regular Blazars would dominate the 
spectrum beyond 20 MeVjj,0j,|(|.) 

If desired, one might try to improve agreement with 
the data in the following way. The continuum (com- 
bined with the other contributions) can fully account 
for the SMM and COMPTEL data, if the clumping fac- 
tor is twice as large as predicted by the NFW profile. 
This could be easily done within uncertainty in our un- 
derstanding of the structure of dark matter halos: for 
example, a slightly steeper profile, or the presence of 
substructure |32j . However, a larger clumping factor also 
increases 511 keV line emission by the same amount, 
which would exceed the HEAO-1 and SMM data. How 
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do we reduce line emission independent of continuum? 
The line emission is suppressed by up to a factor of 4, 
if e + e~ annihilation occurs predominantly via positron- 
ium formation. Once formed, a positronium decays into 
either two 511 keV photons or three continuum photons. 
As the branching ratio of the former process is only 1/4, 
line emission is suppressed by a factor of 4 if all of annihi- 
lation occurs via positronium formation. If a fraction, /, 
of annihilation occurs via positronium, then line is sup- 
pressed by 1 - 3//4 19]; thus, we can cancel the effect 
of doubling the clumping by requiring that 2/3 of line 
emission be produced via positronium. Figure [3] shows 
our "best-fit" model, which assumes (a) nix = 20 MeV, 
(b) the mean clumping factor is twice as large, and (c) 
line emission is solely produced via positronium (/ = 1). 
Note that this is a reasonable extension of the minimal 
model and makes the model more realistic: we know from 
simulations that there must exist substructures in ha- 
los. Some fraction of line emission must be produced via 
positronium, as it has been known that more than 90% 
of 511 keV emission from the Galactic center is actually 
produced via positronium formation [^3, 0] . While the 
model seems to slightly exceed the HEAO-1 and SMM 
data at low energy, we do not take it seriously as the dis- 
crepancy would be smaller than the uncertainty of the 
AGN model. The AGN model presented here assumes a 
high-energy cut-off energy of E cut — 0.5 MeV 3j. Since 
current data of AGNs in such a high energy band are 
fairly limited, uncertainty in E cut is more than a factor 
of 2. Even a slight reduction in E cut would make our 
model fit the low-energy spectrum. 

The best-fit model is consistent with and supported 
by all of the current observational constraints: it fits 
the Galactic 7-ray emission as well as the cosmic 7-ray 
emission. It might also account for a small difference 
between theory and the experimental data of the muon 
and electron anomalous magnetic moment [35|. We stress 
here that, to the best of our knowledge, all of these data 
would remain unexplained otherwise. There is, however, 
one potential conflict with a new analysis of the SPI 
data by |2fj, which shows that a NFW density profile 
does provide a good fit to 511 keV line emission from 
the Galactic center, as opposed to the previous analysis 
by 0], which indicated a shallower profile than NFW. 
This new model would have much higher dark-matter 
clumping and require a substantially (more than an or- 
der of magnitude) smaller annihilation cross-section than 
(av) ~3x 10" 26 cm 3 s^ 1 to fit the Galactic data. Is our 
Galaxy consistent with NFW? This is a rather compli- 
cated issue which is still far from settled (e.g., [13, HH), 
and more studies are required to understand the precise 
shape of density profile of our Galaxy. If our Galaxy 
is described by a steep profile such as NFW, then the 
dark matter annihilation probably makes a negligible 
contribution to the 7-ray background, unless dark matter 
clumping is significantly increased by substructure |32j] . 
compensating a small cross section. On the other hand, 
if it were confirmed that our Galaxy has a shallow density 



profile and the contribution of the dark matter annihila- 
tion to the 7-ray background is negligible, it would be 
difficult to explain the Galactic 7-ray signal solely by an- 
nihilation of light dark matter particles. 

As shown in Figure |3J dark matter annihilation pro- 
duces a distinctive 7-ray spectrum at 0.1-20 MeV. More 
precise determinations of the cosmic 7-ray background 
in this energy band will undoubtedly test our proposal. 
If confirmed, such measurements would shed light on 
the nature of dark matter, and potentially open a win- 
dow to new physics: one implication is that neutralinos 
would be excluded from a candidate list of dark matter. 
Phcnomenologically, our model may be parameterized by 
four free parameters: (1) dark matter mass, m-x, (2) a 
dark matter clumping factor at present, Cx(0), (3) red- 
shift evolution of clumping, /?, and (4) a positronium 
fraction, /. When more precise data are available in the 
future, it might be possible to perform a full likelihood 
analysis and constrain properties of dark matter particles 
as well as dark matter halos. 

Finally, the angular power spectrum of anisotropy of 
the 7-ray background at 1-20 MeV would also offer a 
powerful diagnosis of the detected signal (see [3^] for the 
contribution from Type la supernovae). Our model pre- 
dicts that the angular power spectrum should be given by 
the trispectrum (the Fourier transform of the four-point 
correlation function) of dark matter halos projected on 
the sky, as the signal is proportional to p 2 . More specifi- 
cally, the power spectrum should follow precisely that of 
the dark matter clumping factor. More high-quality data 
of the cosmic 7-ray background in this energy band are 
seriously awaited. 
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Note added in proof - A recent article by Rasera et 
al. (0) argues that our predictions for the 7-ray back- 
ground from the redshifted 511 keV line 

(0) 

were too 

large because annihilation of electrons and positrons can- 
not take place in halos less massive than ~ 1O 7 M0, in 
which baryons cannot collapse. While this effect reduces 
the intensity of the line contribution, it does not affect the 
continuum emission (i.e., the internal bremsstrahlung) , 
as the continuum emission is produced before annihila- 
tion. Since the major contribution to the 7-ray back- 
ground at 1-20 MeV comes from the continuum emission, 
our conclusion in this paper is not affected by the results 
in H3- 
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